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Abstract 

The crossed Andreev reflection in a hybrid nanostructure consisting of a d-wave superconductor and two quantum 
wires is theoretically studied. When the (110) oriented surface of the superconductor is in contact with the wires 
parallel and placed close to each other, the Andreev bound state is formed by the crossed Andreev reflection. The 
conductance has two resonance peaks well below the gap structure in the case of tunnel contacts. These peaks 
originate from the bonding and antibonding Andreev bound states of hole wave functions. 



1. Introduction 

Andreev reflection at the interface of a normal 
metal and a superconductor (SC) [1,2] is one of 
the fundamental consequences of superconductiv- 
ity. This phenomenon corresponds to the scatter- 
ing process of an incoming electron from the nor- 
mal side being reflected as a hole, thereby creating 
a Cooper pair in the condensate. In a tunnel junc- 
tion of a normal metal and a (110) oriented c?-wave 
SC, tiie zero bias conductance peak appears due 
to the formation of the Andreev bound state near 
the interface [3] . 

When a single quantum wire of a single conduct- 
ing channel is in contact with the (110) oriented 
surface of a d-wave SC, the Andreev reflection is 
completely suppressed due to the quantum me- 
chanical diffraction of electron waves at the narrow 
opening [4] . A basic question arises what happens 
if a second quantum wire is placed parallel and 
close to the first one (see Fig. 1). For an electron 
incident from one of the wires into SC, we expect 
that the Andreev hole is reflected back into another 
wire due to the non-local effect called the crossed 
Andreev reflection (CAR) [5,6,7,8,9,10]. Recently, 
crossed Andreev reflection has been observed in an 
s-wave SC with two ferromagnetic wires [11]. 

In this paper, we explore the quantum- 



interference effect due to the crossed Andreev 
reflection in a rf-wave SC with two quantum wires 
parallel and close to each other. It is shown that 
the resonance peak in the conductance is split into 
two sharp peaks at low energies when the bar- 
rier potential of the contact is sufficiently large. 
The lower and higher energy peaks correspond to 
the bonding and the antibonding Andreev bound 
states of quasiparticle holes. 



2. Model and formulation 

We examine the quantum transport in a hybrid 
nanostructure of a c?-wave SC and two quantum 
wires in Fig. 1. A two-dimensional (2D) c?-wave SC 
occupies the right half space, and two quantum 
leads 1 and 2 of width w are connected to SC at 
y = ±L/2. The wave functions of quasiparticles 
(QP) with excitation energy E in the electrodes are 
calculated from the Bogoliubov-dc Gennes equa- 
tion. For simplicity, the Fermi wave number kp and 
the effective mass m are common for all electrodes, 
and the amplitude of the gap function is uniform 
in SC and vanishes in the wires. Our calculation is 
restricted to the case where only the lowest sub- 
band is occupied by electrons or holes in the wires. 
When an electron with energy E and wave num- 




Fig. 1. Schematic diagram of two normal-conducting quan- 
tum wires in contact with a d-wave superconductor. 

ber ki — [2mE + fc| — (vr/w)^]^/^ is incident from 
lead 1 into SC, the wave functions in leads 1 and 
2 are given by '^i{x,y) = (pi{x)x{y — L/2) and 
*2(x, y) = Mx)xiy + L/2) with 
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and x{y) — (2/w)^/^ sin [(tt/w) {y + w/2)], where 
rW and r\\ are the amplitudes of the normal re- 
flection (NR) and the Andreev reflection (AR), re- 
spectively, while r^l and are those of the crossed 
normal reflection (CNR) and the crossed Andreev 
reflection (CAR), respectively. A similar treatment 
is made for an incident electron from lead 2. Since 
i?, A < ep, we put fci « /cf[1 - (tt/Zcfu;)^]^/^ in 
the following. 

To carry out the calculation analytically, 
we make the Andreev approximation following 
Blonder-Tinkham-Klapwijk (BTK) [2] to calcu- 
late the conductance, and put the wave function 
of SC in the form 
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where tf^{py) and t^^iPy) are the transmis- 
sion coefficients, r± — A±/{E + n±), il± — 
^/W^^JA^, and A± = Aocos2(6'Ta) [3,4], 
a being the angle between the (100) axis of SC 



and the normal to the interface (see Fig. 1), and 
6 = si'n~^ {py / kp) is the propagation angle relative 
to the X axis. 

In the following, we focus on the (110) ori- 
ented surface of SC {a — 7r/4) as shown in Fig. 1. 
In this case, A± — ± A with A = Aq sin 26, 
and n± = Q ^ y/E^ - A2. The barrier po- 
tential at the interface between the wires and 
SC is taken into account by the (5-function-type 
potential with amplitude {h^ k-p / 2m) Z , Z be- 
ing a dimcnsionless parameter of barrier height 
[2]. The boundary conditions for the wave func- 
tions at the interfaces are ^'^(0,1/) — 5*^(0,?/) 
and [d^^s{x,y) - d^'i!,{x,y)]^=o = kFZ^'i{0,y) 
{i = 1, 2) appropriate for the S{x) potential. The 
matching technique to the boundary conditions 
[12] yields the reflection coefficients 
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withfci = ki/kF,D± = {ki+T)^-{g,±iZ)^-gl 



(7) 



^= / g§Vi-(p/fcF)V^(p), 

— k-p 

gc^ j ^§V1-(pAf)V'(p)cos(pL), (8) 



dp A 



gs^i J g-Vl-(pAF)V'b)sin(pL), (9) 

where (p{p) — {p\x) is the overlap integral of xiu) 
and e'Py-. 



tp{p) ~ \/^w/tt^cos(j)w/2)/ [1 — (pw/n) 
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Note that t/c and C/^ represent the non-local cou- 
pling and are dependent on L. In the limit of L — > 
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oo, where the two contacts are independent (Qc = 
Gs = O),onehasr^f = {ki-T-iZ)/{ki+T+iZ), 
rfl = r^l = rf^ = 0, indicating the complete sup- 
pression of AR in a single quantum wire with a sin- 
gle transverse mode [4] . The reflection coefficients 
r|| , , 'f'21 , and r'^ for an incident electron from 
lead 2 are obtained from those in Eqs. (3)-(6) by 
the replacement 1 <-> 2 and L — > — i. 



3. Results and discussion 

When the bias voltage V is applied to the two 
leads, the conductance G at zero temperature (T = 
0) is calculated by the formula 
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where h is the Planck constant and E — eV. Fig- 
ure 2 shows the conductance G vs V for kpw — 
4, kpL — 8, and different values of Z. The con- 
ductance is normalized to the normal state value 
Gn ^ G(A = 0). For a low barrier potential of 
small Z, the conductance decreases monotonically 
with decreasing eV below Aq . As the barrier poten- 
tial becomes higher, a peak structure appears well 
below Aq and shifts towards lower eV, developing 
the double peak structure with increasing Z. If the 
conductance is plotted as a function of normalized 
voltage ZV as shown in the inset, the conductance 
peaks fall into the same position, indicating that 
the resonance peak positions are scaled by 1/Z. 

Let us examine the origin of the double peak 
structure in the conductance in the tunneling case 
(Z> 1). Fori; < Ao,J^,^^c, and in Eqs. (7)- (9) 
have the forms: T w ifAo/E, Qc ~ igc^o/E, and 
Gs = igs^o/E, where / = {ET/iAo)E~>o, gc = 
(£'5c/«Ao)£;_o, and gs = {EGs/iAo) are energy 
independent quantities, so that the Andreev re- 
flection coeflacients are calculated in the resonance 
forms 
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Fig. 2. Normalized conductance G/Gm as a function of 
bias voltage V for different values of interfacial barrier 
parameter Z. Inset shows G/Gm vs normalized voltage 
ZV. The circles in the inset are calculated from Eq. (14). 

where g = g,/2/, i?± - (/ ± \gc\ + g^j2f)Ao/Z, 
and 7 = {ki / k-p) f Aq / Z'^ . The other coefficients 
are rff « — 1 — grfj and ~ 0{g^). The reso- 
nance energies E± and intensity g have damped- 
oscillation dependence on L in the period of Fermi 
wave length. We note that the resonance positions 
E± and their separation E+ — E- are scaled by 
l/Z, while the line width 7 is scaled by 1/Z^. 
Therefore the line width of the peaks is much 
smaller than the separation for Z ^ 1, and a 
well-separated two-peak structure is formed as 
shown in the inset of Fig. 2. The conductance G is 
approximated by the sum of two Lorentzians 



G 



77/77 



r,(14) 



Gat (e^-S+)2+72 {eV-E-f+T 

where / = ng"^ f Aq/ ATn and Tn is the normal- 
state value of J-. The simple formula (14) repro- 
duces the numerical result in the inset of Fig. 2, 
if the calculated values (/ = 0.405, gc = —0.036, 
and gs = -0.0345) are used in Eq. (14). Note that 
the peak height of G/Gpf increases in proportion 
to Z\ 

To elucidate the formation of the Andreev bound 
states, we derive an analytical formula of the QP 
wave function ^>s for Z ^ 1 and around E — E± 
in SC: 
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where, except very close to the interface (x ~ 0), 
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with q = y/kl - p2 aj^fj ^ 2Ao(pg/fc|)- At the 
lower resonance energy , is dominated by the 
first term whose electron (hole) wave function ipt. 
{ip'^) is an odd (even) function of y. At the higher 
resonance energy is dominated by the sec- 

ond term whose electron (hole) wave function "0^ 
(ip'l) is an even (odd) function of y. These results 
indicate that the electron (hole) wave functions at 
the lower and higher resonance energies have dif- 
ferent parity with respect to y. 

Figure 3 shows the mapping of the absolute 
squares |V'±(a;,y)P of QP holes at the resonance 
energies for kpw = 5, kpL = 11, and Z = 50. We 
see that the hole wave functions are strongly local- 
ized with large peaks near the contacts due to the 
formation of the Andreev bound states, and that 
the waves are decayed along the (010) and (010) 
directions into which a Cooper pair is created. It 
is noteworthy that the Andreev holes are reflected 
back into leads 1 and 2 in phase a,t E = E- and out 
of phase at E — E+ as seen in Eqs. (12) and (13). 

In summary, we have studied the quantum- 
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Fig. 3. Mapping of Andreev bound states on the xy plane, 
when an electron is incident from lead 1 to a d-wave SC. 
(a) and (b) are the absolute square of hole wave functions, 
li/i^^^l^ and IV'+P, at resonance energies (bonding) and 
E-l- (antibonding) . 



interference effect caused by the crossed Andreev 
reflection in a hybrid nanostructure made up of 
a d-wave SC and two quantum wires. When the 
(110) oriented surface of SC is in contact with 
the wires via tunnel barriers, the Andreev bound 
states are formed at low energies due to the crossed 
Andreev reflection, and two sharp resonance peaks 
appear in the conductance well below the super- 
conducting gap structure. The lower and higher 
conductance peaks correspond to the bonding and 
the antibonding Andreev bound states of the hole 
wave functions. The result predicted here is pos- 
sibly observable in a d-wave SC with conducting 
molecular wires, self-organized nanowires, or a 
STM tip with two atomic point contacts. 
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